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ORIGINAL ARTICLE 

Quantification of subclonal distributions of recurrent genomic 
aberrations in paired pre-treatment and relapse samples from 
patients with B-cell chronic lymphocytic leukemia 

SJL Knight 1 ' 2 , C Yau 3 , R Clifford 4 , AT Timbs 4 , E Sadighi Akha 2 , HM Dreau 4 , A Burns 1 ' 4 , C Ciria 5 , DG Oscier 6 , AR Pettitt 7 , S Dutton 5 , 
CC Holmes 2 ' 3 , J Taylor 1 ' 2 , J-B Cazier 2 and A Schuh 1 ' 4 

Genome-wide array approaches and sequencing analyses are powerful tools for identifying genetic aberrations in cancers, 
including leukemias and lymphomas. However, the clinical and biological significance of such aberrations and their subclonal 
distribution are poorly understood. Here, we present the first genome-wide array based study of pre-treatment and relapse 
samples from patients with B-cell chronic lymphocytic leukemia (B-CLL) that uses the computational statistical tool OncoSNP. 
We show that quantification of the proportion of copy number alterations (CNAs) and copy neutral loss of heterozygosity 
regions (cnLOHs) in each sample is feasible. Furthermore, we (i) reveal complex changes in the subclonal architecture of 
paired samples at relapse compared with pre-treatment, (ii) provide evidence supporting an association between increased 
genomic complexity and poor clinical outcome (iii) report previously undefined, recurrent CNA/cnLOH regions that expand 
or newly occur at relapse and therefore might harbor candidate driver genes of relapse and/or chemotherapy resistance. 
Our findings are likely to impact on future therapeutic strategies aimed towards selecting effective and individually tailored 
targeted therapies. 
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INTRODUCTION 

Genome-wide array approaches and sequencing analyses are 
powerful tools for identifying genetic aberrations in cancers, 
including leukemias and lymphomas. 1 10 B-cell chronic lympho- 
cytic leukemia (B-CLL) is the most common adult leukemia in the 
Western World. 11,12 It is characterized by a chronic relapsing 
course and the development of chemotherapy refractoriness, 
leading to death in a significant subset of patients. There is, 
therefore, a clinical need to develop novel treatment strategies 
that overcome chemotherapy resistance. 

B-CLL can be subclassified according to immunoglobulin heavy 
chain variable gene sequence (IGHV) homology into 'mutated' or 
'unmutated' cases, reflecting good or poor prognosis, respectively. 
B-CLL shares many of the molecular characteristics seen in other 
forms of cancer and recurrent copy number alteration (CNAs) with 
relevance to prognosis have been described, including trisomy 
of chromosome 12 (16%) and loss of chromosomal regions 
17p13.1 {TP53; 7%), 13 " 17 1 1q22.3 {ATM; 18%) 18 " 19 and 13q14.3 (DLEU2 
and/or miR-15a/16-1; 55%). 20 " 30 In addition, single-nucleotide poly- 
morphism (SNP) and comparative genomic hybridization arrays 
have allowed novel recurrent genomic abnormalities to be 
identified 31 " 33 that have been linked to prognosis in B-CLL. 30 ' 34 " 38 

Current models of cancer progression are based on the concept 
that tumors are subject to the Darwinian process of evolution and 
selection 39,40 Recent studies in acute lymphoblastic leukemia have 



provided pivotal insights into the complex sequence of events 
during leukemogenesis, showing that the initiating mutation is 
followed by CNAs that drive the emergence of frank leukemia. 39 " 42 
Together, these data imply that at least some CNAs/copy neutral loss 
of heterozygosity regions (cnLOHs) are likely to be involved in driving 
leukemia progression and therefore might contribute to relapse. 

This led us to hypothesize that subclones containing driver 
CNAs/cnLOHs would newly occur or expand in relapse samples 
compared with samples taken before treatment and would be 
recurrent within our patient cohort. In contrast, random passenger 
mutations would remain unchanged or decrease/disappear in 
paired pre-treatment/relapse samples and would not be observed 
recurrently in different patients within the cohort. If correct, then 
we would anticipate that the identification of driver genes within 
recurrent or emerging/expanding regions of CNA/cnLOH might, in 
the longer term, have the potential to inform the design of novel 
therapies aimed at treating relapsed B-CLL. 

In the present study we tested our hypothesis by systematically 
tracking the presence and subclonal distribution of CNAs/cnLOHs 
in patients before treatment and at subsequent relapse. To 
achieve this, we used high-resolution SNP array technology. We 
chose the newly developed computational statistical tool 
OncoSNP 10 that provides quantitative measures of cell admixture 
on a per-SNP level. OncoSNP was selected in preference to other 
dedicated cancer tools for SNP array analysis, 43 " 47 because our 
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data were generated from samples at two time points only and 
our study was specific to B-CLL, where spatial heterogeneity (for 
example, samples biopsied from different parts of a tumor, 
metastases, etc) is not applicable (see Supplementary Information, 
online). Furthermore, we were interested not only in the 
accumulation of genomic aberrations over time, but also 
particularly in the expansion of distinct subclones. Thus, OncoSNP 
best suited our study where (i) samples differed from solid tumors 
in being generally diploid with comparatively fewer CNAs (ii) 
normal cell contamination could be minimized and (iii) we had a 
specific interest in subclonal evolution. 

PATIENTS AND METHODS 

For a detailed account of the Materials, Patients and Methods used 
in this study please refer to the Supplementary Information online. 

Patients and samples 

All patients gave written informed consent in accordance with the 
Declaration of Helsinki. In all, 135 samples from 93 patients were 
included in our cohort. For 42 of the 93 patients, paired pre- 
treatment and relapse samples were available. Only samples from 
patients with lymphocytes contributing >90% of the total white 
blood cell count were included in the analysis. DNA was extracted 
from vital frozen cells using the QIAamp DNA Midi Kit (Qiagen, 
Hilden, Germany) according to the manufacturer's protocol. 

Fluorescence in situ hybridization (FISH) 

Interphase cells from CLL cultures were analyzed by FISH using the 
Vysis CLL FISH panel probe set (Abbott Molecular, Illinois, IL, USA) 
according to manufacturer's instructions. 

Arrays 

Hybridization to lllumina Genome-wide SNP chips was performed 
according to manufacturer's protocols found on registration at 
http://www.illumina.com/products/human1m_duo_dna_analysis_ 
beadchip_kits.ilmn and at http://www.illumina.com/support/array/ 
array_kits/humanomni1-quad_beadchip_kit/documentation.ilmn. 
The data were processed using GenomeStudioV2009.2 (lllumina, Inc., 
San Diego, CA, USA) and then analyzed using OncoSNP v1.0 (see 
below). For visual comparisons, the data were processed also using 
Nexus 5 Discovery Edition (BioDiscovery, Inc., El Segundo, CA, USA). 

OncoSNP analysis 

Detailed methods have been described previously 10 and further 
details are given in the Supplementary Information online. To help 
identify and exclude germline CNVs of unlikely relevance from our 
data, we excluded SNP and monomorphic copy number probes in 
known germline CNV regions from the OncoSNP analysis and 
we also made use of the database of genomic variants 48 and 
the Wellcome Trust Case Control Consortium cohort data (see 
www.wtccc.org.uk). 

TP53 mutation analysis 

TP53 was screened for mutations using four high-resolution 
melting assays covering exons 5-8. Following high-resolution 
melting, results were analyzed on the high-resolution melting 
module of the Corbett Rotor-Gene 6000 software 1.7 (available for 
download from QIAGEN, http://www.qiagen.com/corbett/support/ 
default.aspx). Positive PCR products were purified and then 
sequenced from both strands on the CEQ 8000 (Beckman Coulter, 
Fullerton, CA, USA). 

IGHV mutation analysis 

To identify clonal rearrangements of the IGHV gene and determine 
the somatic mutation status of the variable (V) gene sequence in 
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patients we used the IGH Somatic Hypermutation Assay v2.0 
(Invivoscribe, La Coutat, France) according to the manufacturer's 
instructions. 

RESULTS 

A flowchart outlining the patient numbers and overall study 
design is given in Supplementary Figure 1 online. 

Our preliminary experiments demonstrated the importance of 
selecting an array platform capable of detecting CNAs and cnLOH 
(see Supplementary Information, Supplementary Table 1 and 
Figure 2 online) as we found regions of cnLOH signposting genes 
that might carry mutations relevant for treatment choice and that 
would otherwise have been missed using a non-SNP-based 
platform (for example, TP53, see Supplementary Information, 
Supplementary Table 2 and Figure 3 online). We also demon- 
strated in these experiments that OncoSNP provided accurate 
quantification of CNAs, based on the strong correlation between 
OncoSNP results and conventional FISH analysis (Supplementary 
Table 3 and Figure 4 online). 

Genome-wide analysis of pre-treatment samples 
For our main study, we focused initially on the analysis of pre- 
treatment samples from 93 patients undergoing regular follow-up 
at our institutions. This cohort was representative of CLL cohorts 
with a higher proportion of unmutated IGHV genes and/or 
treatment resistant cases that are referred to tertiary centers for 
treatment consideration (Table 1). A total of 80 patients had one 
or more CNAs/cnLOHs in well-established regions of prognostic 
importance - 17p13.1, trisomy 12, 13q14.3, 11q22.3 (Table 1, 
Supplementary Table 2, Figures 2, 3 and 5 online). CNAs involving 
chromosome arms 8p, 9p and 10q developed exclusively in 
patients with unmutated IGHV. 

We identified 58 previously unreported recurrent CNAs/cnLOHs 
(Table 2) (see Supplementary Information online for recurrence 
criteria). We established the minimally deleted regions (MDRs)/ 
minimally overlapping regions (MORs) of the recurrent CNAs/ 
cnLOHs and revealed interesting candidate genes including those 
implicated in B-cell maturation (for example, BLIMP1, NFkB2, TLR4, 
and CREBBP), tumor progression (for example, RND3, RHOT1; 
RHDBL; RAB20, and TRAP1), DNA damage response (for example, 
RIFI, JP53, AJG4D, and AJG5), tumor suppression (Tor example, 
NMI, CNOT7, PDGFRL, FGF20, PI3K, and FOG2) and familial CLL (for 
example, SP140L, SP100) (Table 2). Figure 1 shows examples of four 
of these novel recurrent regions of interest that span only a few 
candidate genes, potentially facilitating the identification of driver 
mutations. Sequencing of the TP53 gene, the known candidate 
underlying poor prognosis on 17p13.1, revealed mutations in 14 
out of 93 patients and deletions in 1 1. However, exon sequencing 
of two of the emergent candidates, BLIMP1 and ATG5 revealed no 
mutations in 40 patients. 

In order to establish the clinical significance of our findings, we 
correlated our array results with measures of clinical prognosis 
and risk, including IGHV mutation status and clinical risk scores 
defined by necessity for treatment intervention, progression free 
survival and chemotherapy resistance (see Supplementary Infor- 
mation, online). Genomic complexity was defined either by the 
presence of three or more CNAs ^ 20 kb and/or cnLOHs ^ 2 Mb in 
addition to, or other than, the known CNA regions assayed 
routinely by FISH, or by a total length of CNAs ^5 Mb. Within the 
93 pre-treatment samples, patients with unmutated IGHV genes 
showed a statistically significant higher total number of CNAs/ 
cnLOHs (50.0% ^3 CNAs/cnLOHs) than patients with mutated 
IGHV genes (16.0% ^3 CNAs/cnLOHs) (Cochran -Armitage trend 
exact test, P-value = 0.0032, Supplementary Table 4 online) and 
patients with del17p/7P53 mutation showed a higher total length 
of CNAs/cnLOHs (90.9% >5 Mb CNAs/cnLOHs) (Cochran -Armitage 
trend exact test, P-value = 0.0081, Supplementary Table 5, online). 
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Table 1. Characteristics of patients included in the study 


Pre-treatment analysis (n = 93) 


Patient Characteristics 




Patient numbers 


At time of blood sample: 






Mean age (range) 




66.75 (36-87) 






years 


Females 




28 


Males 




65 


dell 1q22.3 




17 


dell 3q 14.3 




47 (28 as sole 






abnormality) 


Trisomy 12 




14 


del17p13.1 




8 


None of these 




26 


TP53 bi-allelic mutations 




3 


TP53 mono-allelic mutation 




3 


TP53 mutation/dell 7p 




8 


TP53 wt/dell7p 




1 


TP 53 mutation negative 




68 


TP53 mutation not done 




10 


Mutated IGHV genes 




25 


Unmutated IGHV genes 




42 


Not known 




26 


W and W (CRS = 0) 




16 


PFS>1 year (CRS = 1) 




22 


PFS<1 year (CRS = 2) 




20 


Refractory (CRS = 3) 




18 


FU<1 year (CRS = N/A) 




15 


NK (CRS = N/A) 




2 


Treatments following blood sampling 






Chlorambucil 




16 


Chlorambucil + Rituximab or 




3 


Chlorambucil + Ofatumumab 






Bendamustine single agent 




1 


Purine analogue combination 




1 2 


Purine analogue combination + 




11 


Rituximab 






Ofatumumab single agent 




1 


Alemtuzumab 




1 


Chemotherapy refractory 




18 


NK 




14 


Number of CNAs/cnLOHs 






None 




13 


<3 




39 


^3 




41 


Length of CNAs /cnLOHs 






<1Mb 




32 


1 -5Mb 




1 1 


>5Mb 




50 


Paired (Pre-treatment and Relapse) analysis (n 


= 42) 


Patient Characteristics 


Patient numbers 


Clonal evolution 






(II — ZD J 


dell1q22.3 


5 


4 


del13q14.3 


21 (13 as sole 


7 




abnormality) 




Trisomy 12 


7 


2 


del17p1 3.1 /TP53 mutation 


5 


5 


None of these 


12 


5 


Mutated IGHV genes 


13 


7 a 


Unmutated IGHV genes 


25 


14 b 


NK 


4 


4 


PFS>1 year (CRS = 1) 


19 


10 


PFS<1 year (CRS = 2) 


18 


10 


Chemotherapy refractory (CRS = 3) 


5 


5 


Abbreviations: CNA, copy number alteration; cnLOH, copy neutral loss 


of heterozygosity; CRS, clinical risk score; FU, follow-up after treatment; 


IGVH, immunoglobulin variable heavy chain gene; NK, not known; PFS, 


progression free survival; W and W, Watch and Wait. a 6/7 had % change 


only. 11/14 had newly emerging CNAs/cnLOHs. 





In our small cohort, large Type II 13q deletions were not associated 
with a worse clinical outcome (Supplementary Tables 6 and 7 online). 

Furthermore, patients with high clinical risk scores had a higher 
number of CNAs/cnLOHs than patients with low clinical risk scores 
(Kruskal-Wallis test, P-va I ue = 0.0016) and showed a greater 
overall length of CNAs/cnLOHs than patients with low clinical risk 
scores (Kruskal-Wallis test, P-value = 0.0002) (Supplementary 
Table 8, online). Importantly, 16 patients had ^3 CNAs/cnLOH 
but no 17p13.1 or 11q22.3 loss and would not have been picked 
up by FISH as a poor risk group. 



Genome-wide comparison of paired pre-treatment/relapse 
samples 

Next, we focused on the analysis of the 42/93 patients for whom 
paired samples (that is, both pre-treatment and relapse samples) 
were available (Table 3). We hypothesized that in addition to 
recurrent CNAs, newly occurring or expanding genomic aberra- 
tions would also represent potential drivers of disease progres- 
sion. The percentages of subclones carrying CNAs/cnLOH regions 
were calculated from the SNP data using OncoSNP. 10 Only one 
case with mutated IGHV genes (CLL086) had additional CNAs at 
relapse. The other changes were percentage changes of pre- 
existing CNAs. By contrast, 1 1 unmutated IGHV genes cases had 
complex evolution with emergence of one or more additional 
CNAs. When we looked at total length of CNAs/cnLOHs, patients 
without clonal evolution showed a statistically significant lower 
total length of CNAs/cnLOHs (58.8% <1 Mb) than patients with 
clonal evolution (68.0% >5 Mb)(Cochran-Armitage trend exact 
test, P-value = 0.0075; Supplementary Table 9, online). 

All 13q14.3 deletions were present before treatment and of the 
21 patients manifesting these; seven exhibited ^10% changed 
proportions of subclones carrying the deletion at relapse. In two 
cases, (CLL081 and CLL107), the proportion of subclones with the 
del 13q 14.3 appeared decreased at relapse and both exhibited 
additional CNAs elsewhere in the genome. For another two cases 
(CLL080 and CLL096), with deletions > 1 Mb, the proportion of 
subclones with the 13q14.3 deletion had increased over time; 
both had isolated chromosome 13 anomalies (Figure 2). Further- 
more, in CLL071 and CLL080, there was extension of the 1 3q1 3.4 
deletion in one allele (Supplementary Figure 6 online, Supple- 
mentary Table 5). Overall, the results suggest that extension of the 
deleted 13q14.3 locus may confer a clonal survival advantage in 
the minority of cases without additional CNAs/cnLOHs. By 
contrast, for patients with complex genomic aberrations, clones 
containing dell 3q1 4.3 are outcompeted. 

Patients with del17p13.1/TP53 mutations were most likely to 
acquire additional CNAs at relapse (5/5). In addition, three patients 
without dell 7p1 3.1 at diagnosis developed del17p and other 
CNAs at relapse (CLL084 and CLL108, CLL145). Patient CLL081 was 
a clear example of this, manifesting a 17p13.1 deletion and 
multiple other CNAs pre-treatment. At the time of relapse, five 
CNAs had altered in proportion (Figure 3) and two events 
(del2q33.1-q37.1 and dell 6p1 3) were newly identified. Further- 
more, a number of the CNAs identified encompassed newly 
identified MDRs/MORs (Table 1) and included the genes LDB-1 and 
NFkB2 (Wq2432 MDR), SP140L and SP100 (2q37.1 MOR) and 
BTBD12, DNASE1, TRAP! and CREBBP f16p13.3 MDR). Finally, the 
OncoSNP feature of defining proportions of subclones carrying 
CNAs/cnLOHs allowed us to use data from matched pre-treatment 
and relapse samples to infer likely subclonal populations and to 
map their evolutionary relationships. Schematic representations of 
the likely clonal architecture at different time-points for patients 
CLL081 and CLL106 are shown in Figure 4 with a further example 
for CLL092 in Supplementary Figure 7 online. These illustrate 
clearly a non-linear, branching subclonal hierarchy in B-CLL with 
multiple ancestral subclones, already present in pre-treatment 
samples and contributing to relapse. In conclusion, using the 
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Table 2. Newly defined regions showing recurrent/expanding CNAs/cnLOHs 


Recurrent CNAs/cnLOH 


No. of 


No. of 


No. of 


MOR 


MOR 


Genes/microRNAs 




pretreotment 


new 


events with 


start (bp) 


end (bp) 






samples 


events 


increased 










showing 


observed 


percentage 










overlap 


at relapse 


at relapse 








loss/gain/cnLOH 1 p36.31 p36.23 


3 


— 


— 


6508255 


7 966180 


Many including THAP3, CAMTA1, 














TNFRSF9/CD137 


loss/cnLOH 2q31.1 


1 


1 


— 


174294144 


174 955141 


SP3, OLA1/GTPBP9, SP9, CIR1/CBF1 


loss/cnLOH 2q36.2 


3 


1 


— 


225 091 358 


225 626330 


CUL3, DOCK10 


loss/cnLOH 2q37.1 


4 


1 


— 


230940929 


231 244281 


SP140L, SPW0 


loss/gain 3p25.3 


2 a 


— 


— 


9610471 


10363 762 


Many including CIDEC 


loss/gain 3p24.3 


2 a 


— 


1 


21 592 730 


22 039016 


ZNF385D 


loss/gain 3p21.31 


2 a 


— 


— 


46827 053 


47426 525 


Many including PTPN23 


loss/gain 3p21.1 


2 a 


1 


— 


52 038099 


52 553 933 


STAB 1 /CLEVER 1/FEEL 1, NT5DC2, C3orf78 


loss 4q13.1-q13.2 


2 


— 


— 


64455 242 


66975 683 


TECRL, LOC401134, EPHA5, hsa-mir-1 269, 














LOO 001 44602 


loss 4q13.2q13.3 


2 


— 


— 


69847111 


71 783 789 


Many including IGJ, UTP3 


loss 5q21.3q22.1 


3 


— 


— 


106355 299 


110412993 


EFNA5, FBXL17, FER, PJA2, MAN2A1, 














LOd 00289673, TMEM232, SLC25A46 


loss 5q22.2q22.3 


3 


— 


— 


112371 443 


113 463403 


DCP1, MCC, TSSK1B, YTHDC2 


loss/gain 6p25.3-p25.2 


2 


2 


— 


3138439 


4 294807 


TUBB2B, PSMG4, SLC22A23, C6orf145, 














FAM50B, PRPF4B, PEO 


loss 6q16.3 


4 


2 


— 


101 257101 


101 445 547 


ASCC3 


loss 6q21 


4 


— 


— 


106208162 


106 917 920 


BLIMP-1, ATG5 


loss 6q21 


3 


1 


— 


109405 342 


109444628 


SESN1 


loss7q36 


2 


1 


1 


156694855 


156 824049 


UBE3C, KIAA0010, MRJ, DNAJB6 


loss 8p22p21.2 


3 


2 


— 


16324 942 


20039 997 


Many including CNOT7, FGF20, MTUS1, 














PDGFRL 


loss 8p12 


2 


2 


— 


33488001 


36433 467 


c8orf41, RNF122, DUSP26, UNC5D 


loss 8p11.22p11.21 


2 


1 


- 


39 796 895 


40860 746 


ADAM2, ID01, ID02, C8orf4, ZMAT4 


loss 8p11.2p11.1 


2 


1 


— 


40947 586 


43 533112 


SFRP1, GOLGA7, AGPAT6, NKX6-3, ANK1, 














MIR486, MYST3, AP3M2, PLAT, IKBKB, VDAC3, 














SLC20A2, C8orf40 CHRNB3, CHRNA6, THAP1, 














HOOK3, RNF170, FNTA, SGK196, POTEA 


loss/gain/cnLOH 8q23.1 


3 


1 


1 


106323 747 


106 959523 


ZFPM2 (FOG2) 


loss/gain 9p22.3-p24.1 


5 


- 


1 


11 986163 


14313346 


TYRP1, MPDZ, NFIB 


loss/cnLOH 9q33.1 


3 


— 


— 


119325172 


121 100178 


TLR4, DBC1 


loss/cnLOH 9q33.1-q33.2 


3 


— 


— 


121 461 281 


122 566552 


MIR147, CDK5RAP2, MEGF9, FBXW2 


loss/cnLOH 9q34.13-q34.3 


3 


— 


— 


134588424 


139101 448 


Many including hsa-mir-126 


loss 10p12.33 


2 


— 


— 


17 900 269 


18398113 


FAM23A, MRC1L1, MRC1, MIR51 1-1, SLC39A12 


loss/cnLOH 10q23.2 


2 


— 


1 


89310000 


89454 971 


PAPSS2/ATPSK2 


loss 10q23.33 


- 


2 


- 


95 437 455 


95 462 005 


FRA10AC1 


loss 10q24.32 


4 


3 


— 


103 708264 


104 851 056 


Many including LDB1; NFkB2; SUFU; TRIM8 


loss /cnLOH 11q12.2q13.1 


2 


— 


- 


61 231 809 


64 790450 


Many including MTA2/MTA 1L 1/PID 


loss/cnLOH11q24.2q24.3 


4 


1 


— 


126353410 


126 814884 


KIRREL3, MIR3167 


loss/cnLOH11q24.3 


4 


1 


— 


127 852 266 


128613424 


ETS1, FLU, KCNJ1, c11orf45, TP53AIP1, 














ARHGAP32 


loss/cnLOH11q24.3 


4 


1 


— 


128805 080 


128 951 253 


BARX2 


loss/gain/cnLOH12p13.33 


2 a 


— 


— 


977 044 


1 959491 


ERC 1/ELKS/RAB6IP2, LOC1 00292680, FBXL14, 


loss/gain/cnLOH1 2p1 3.2p1 3.1 












WNT5B, ADIPOR2, CACNA2D4, LRTM2 


2 a 


— 


— 


10176350 


14248 326 


Many including ETV6 


gain/cnLOH12q13.11q13.12 


4 a 


— 


— 


47 384 207 


48029106 


Many including MLL2 


loss/gain 12q14.3 


2 a 


— 


— 


64809667 


64853 817 


LLPH, S1R, TMBIM4 


loss/gain/cnLOH 12q23.1 


4 a 


— 


— 


99 1 75 388 


99224488 


DEPDC4, SCYL2/CVAK1 04 


loss/gain/cnLOH 12q24.11 


4 a 


— 


— 


109411 836 


109658160 


C12orf24, VPS29, RAD9B, PPTC7,(TAPP2C) 














TCTN1, HVCN1, PPP1CC 


gain/cnLOH12q24.31 


4 a 


— 


— 


119578493 


121 172 890 


Many including ORAI1 and BCL7A 


loss/cnLOH 13q14.11 


4 


— 


1 


39 760000 


39910000 


LOC646982/TTL-B1 


loss/cnLOH 13q14.11 


6 


— 


2 


40389170 


40 761 341 


SUGT1L1, ELF1, WBP4, KBTBD6, MTRF1 


loss/cnLOH 13q14.11 


6 




2 


41 687 912 


41 826016 


DGKH, AKAP1 1/AKAP220 


loss/gain/cnLOH 13q14.2 


32 


— 


1 


47 781 909 


47 832 800 


RBI 


loss/gain/cnLOH 13q14.2 


22 


— 


1 


47 387 659 


47 522 891 


SUCLA2, NUDT 


loss/cnLOH 13q33.3-q34 


5 


1 


— 


108839643 


112014298 


Many including RAB20, ING1, ARHGEF7 


loss/gain 1 5q1 5.1 


2 


- 


- 


38967 655 


40822 359 


Many including NDUFAF1 


loss 16p13.3 


1 


1 




3 1 96 980 


4492425 


Many including ADCY1, TFAP4, GLIS, 














DNAJA3, HMOX2, CDIP 


loss/gain 1 7q1 1 .2 


2 


1 




24979666 


25 506441 


SSH2, EFCAB5, CCDC55, MIR423, MIR3184 


loss/gain 17q1 1.2 


2 


1 




27442 088 


27 692 902 


RHOT1, ARGFXP2, RHBDL3, C17orf75 


loss/gain 17q12 


3 






30 750407 


31 140027 


Many including AP2B1, RASL10B, MMP28 


gain 17q12 


3 






34817 937 


34948 998 


CDK12, CRKR5, MED 1/TRAP220/PBP/PPARBP 


loss/cnLOH 18p11.32p31 


4 


1 


1 


1 


6 572104 


Many including USP14, THOC1 /nuclear 














matrix protein p84/HPR1) 
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Table 2 {Continued) 


Recurrent CNAs/cnLOH 


No. of 


No. of 


No. of 


MOR 


MOR 


Genes/microRNAs 




pretreatment 


new 


events with 


start (bp) 


end (bp) 






samples 


events 


increased 










showing 


observed 


percentage 










overlap 


at relapse 


at relapse 








gain 19p13.2 


2 a 


1 


— 


10071 486 


10926409 


Many including ATG4D, hsa-mir-1 181, 














hsa-mir-1238 


loss 19p13.11 


1 a 


1 


— 


18336675 


19169002 


Many including SSBP4, FKBP8, UBA52, 














CRTC1, BMP family (GDF1, LASS1), ELL 


loss/ LOH loss19q13.12q13.13 


2 a 


— 


— 


41 440 201 


43149438 


Many including many ZNF genes 


gain/cnLOH 20q13.13 


4 


— 


1 


47 292 495 


47 500 772 


DDX27, ZNFX1, SNORD12C, SNORD12B, 














SNORD12, C20orfl99, NCRNA00275, 














AK055386 


loss/cnLOH 20q11.23 


4 




1 


34645 412 


34974 353 


TGIF2, C20orf24, SLA2, NDRG3, DSN1, 














C20orfll7, C20orfll8, SAMHD 7 


loss 20q 11.23 


5 




1 


35 010463 


35 020688 


SAMHD1 


loss/gain 22q12.1 


2 


1 




26427 496 


26490006 


MA/ 7 


loss/gain 22q12.1q12.2 


2 






27 573 410 


27918140 


ZNRF3, C22orf3h KREMEN1 


loss/gain/cnLOH Xq27.3 


3 


1 




142 536091 


144254333 


SLITRK4, SPANXN2, SPANXN1 


Abbreviation: CNA, copy number alterations; cnLOH, copy neutral loss of heterozygosity regions; MOR, minimally overlapping region. All co-ordinates are 


hg36. Shaded blue results show the only genes/microRNAs mapping within the identified MORs. Excluding full trisomies. 
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Figure 1. Four recurrent regions of interest observed in B-CLL patient samples. The boundaries of the MDRs/MORs, 2q37.1, 6q21, 9p22.3p24.1 
and 10q24.32, are indicated by the vertical lines. Red lines show called losses, green lines show gains, purple lines show allelic imbalance and 
mustard lines show cnLOH. In each case the MDR/MOR was verified by visual inspection of Log ratio and B-allele frequency plots. P = Pre- 
treatment, R = Relapse. 



combination of high-resolution SNP array and OncoSNP, we were 
able to detect clonal evolution in 60% of treated patients with B-CLL 
Interestingly, when comparing the regions affected by clonal 
expansion with recurrent regions we noticed a considerable 
overlap (Table 3). Many regions that were both recurrent and 



expanded or newly occurred in relapse samples included a very 
limited number of genes (Table 2, blue shading). These regions 
include genes that have a role in familial CLL {SP140L, SP100), in 
B-cell development and autophagy {BLIMP1 and ATG5, respectively) 
and FOG2, a regulator of phosphatidylinositol 3-kinase, which is 
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Table 3. 


Comparison between CNA and cnLOH events noted in sequential pre-treatment versus relapse samples 3 


Patient 


Known 


IGHV gene 


CR5 


CNAs and cnLOH regions in pre-treatment sample 


Differences at relapse 


no. 


CNAs 


status 








CLL010 


None 


Un mutated 


1 


None noted 


None noted 


CLL023 


None 


Un mutated 


1 


46,XX, arr 4q24(1 07,1 62,437-1 07,237,623)x3 b , 


46,XX, arr 1 0pl 2.23(1 7,833,779-1 8,398,1 13)x1 no 










40% 10p1 2.23(1 7,833,779-1 8,398,1 13)x1, 40% 


longer detectable, -10% 10q24.32q25.2 










1 0q24.32q25.2 (1 03,562,91 5-1 1 2,637,1 85)x1 


(1 03,562,91 5-1 1 2,637,1 85)x1 


CLL035 


Tri12 


Mutated 


1 


47,XX, arr 3p24.3(21,592,730-22,039,016)x1 b , 


None noted 










12(0-1 32,1 46,663)x3 




CLL050 


del 13q 14.3 


Un mutated 


1 


46, XX, arr 6q14.1 (83,1 65,91 7-83,265,032)x1 b , 


None noted 










7p21 .1 (1 8,576,853-1 8,602,061 )x3 b , 1 1 q 1 2.3 












(62,082,265-62,330,255)x1 b , 1 3q 1 4.2q 1 4.3 












(48,608,987-5 1,326,824)x1 




CLL068 


del 13q 14.3 


Mutated 


1 


46,XY, arr 90% 2q22.2q24.1 (144,240,1 56- 


46,XY, arr >90% 2q22.2q24.1 (144,240,1 56- 










1 57,961 ,1 29)x1 , 1 3q 1 4.2q 1 4.3(47,595,821 - 


1 57,961,1 29)x1 










50,072,31 2)x1 




CLL080 


del 13q 14.3 


Mutated 


1 


46,XX, arr 3q25.31 (156,553,1 59-1 56,572,680)x1 b , 


46,XX, arr 80% 1 3q1 4.1 Iql 4.3(40,691, 755- 










30% 13q14.1 1q14.3(40,674,362-49,342,922)x1 b , 


50,580,804)x1 (extended across region) 










- 50% 1 3q 1 4.3(49,342,923-50,332,1 92)x1 b 




CLL083 


None 


Un mutated 


1 


46,XX, arr 4p1 4(39,001, 471 -39,049,059)x3 b , 


46,XY, arr 1 1q25(1 31,1 23,91 0-1 31,1 86,533)x3 b no 










1 1 q25(1 31,1 23,91 0-1 31 ,1 86,533)x3 b , < 1 0% 


longer detectable. 1 2q 14.3 (64,809,667- 










1 2q1 4.3(64,809,667-64,853,81 7)x1 b , 


64,853,81 7)x1 b no longer detectable, 










1 2q23.1 (99,1 75,388-99,224,488)x1 b ,1 0% 


12q23.1(99,175,388-99,224,488)x1 b no longer 










1 6p1 2.2(20,71 0,305-20,790,857)x3 b , 


detectable, 1 6p1 2.2(20,71 0,305-20,790,857)x3 b 










1 7q1 2(34,81 7,937-34,948,998)x3 b 


no longer detectable, 17q 12(34,81 7,937- 












34,948,998)x3 b no longer detectable 


CLL085 


del 13q 14.3 


Un mutated 


1 


46,XY, arr 10q23.1q23.3 1(85,670,734- 


46,XY, arr 10q23.33(95,437,455-95,462,005)x1 b 










89,454,971)x2 hmz, 50% 1 3q1 4.1 1 q21 .31 












(43,934,565-61,368,1 81 )x1 




CLL089 


None 


Un mutated 


1 


a a vv -mtv 1 «ti i/i/i/i qqc i on 1/1/1 a 1 /i q/i 

4o, ay, arr I qz \ .\(\ 44,dd:>, i dU- \ 44,o 1 4,y4z>jxo , 


None noted 










1 q25.2(1 74,495,962-1 74,547,429)x3 b 




CLL090 


tri12 


Mutated 


1 


47,XX, arr 90% 7q33q34(1 34,397,353- 


47,XX, arr >90% 7q33q34(1 34,397,353- 




del 13q 14.3 






139,656,577)x1, 90% 7q36 (156,694,855- 


139,656,577)x1, >90% 7q36(1 56,694,855- 










1 56,824,049)x1 b , 50% 1 1 q1 2.1 qter(60,468,463- 


156,824,049)x1 b , 90% 1 3q 14.3(49,296,907- 










134,445,626)x2 hmz, 80% 12(0-1 32,1 46,663)x3, 


50,580,804)x1 










80% 1 3q 1 4.3(49,296,907-50,580,804)x1 




CLL091 


None 


Mutated 


1 


46,XY, arr 14q24.3q31.1 (77,490,883-80,987,659)x2 
hmz c 


None noted 


CLL093 


del 13q 14.3 


Mutated 


1 


46,XY, arr 80% 6p25.3p25.2-(94,609-4,294,807)x3, 


None noted 










6q 1 4.1 -qter(83,755,489-1 70,753,209)x1 , 70% 












1 0q25.1 (1 09,682,821 -1 09,770,21 5)x3 b , 1 3q1 4.2- 












q 1 4.3(47,556,552-50,61 0,384)x1 , 70% 












20q1 3.1 3(47,292,495-47,500,772)x3 b 




CLL096 


del 13q 14.3 


Mutated 


1 


46,XX arr 40% 13(0-1 14,1 21, 252)x2 hmz, 


46,XX, arr 90% 13(0-1 14,1 21, 252)x2 hmz. 13q14.3 










1 3q1 4.3(49,425,974-49,522,1 41 )x1 , 


(49,425,974-49,522,1 41 )x1 no longer detectable 










1 3q 1 4.3(49,522,1 41 -50,380,71 3)x0 




CLL097 


del 13q 14.3 


Mutated 


1 


46, XY, arr >90% 1 3q1 4.1 1q1 4.3(40,332,062- 


46,XY, arr 90% 1 3q1 4.1 1q1 4.3(40,332,062- 










51,453,289)x1 


51,453,289)x1 


CLL098 


None 


Un mutated 


1 


None noted 


None noted 


CLL099 


None 


Un mutated 


1 


46,XY, arr 70% 1q21 .1(1 44,099,494- 


None noted 










144,757,265)x3 b 




CLL101 


de11q22.3 


Un mutated 


1 


46,XY, arr 1 1q14.1qter(77,080,067-134,445,626)x1, 


None noted 




del 13q 14.3 






60-70% 1 3q 1 4.2q21 .32(47,781 ,909-64,41 0,990)x1 b , 












50% 1 3q21 .32q33.2(64,350,990-1 03,707,225)x3 b , 












70% 22q1 1 .21 qter(1 7,403,754-49,582,267)x3 




CLL107 


de11q22.3; 


Un mutated 


1 


46,XY, arr 50% 3p21. 31 (46,827,053-50,1 33,778)x1, 


46,XY, arr 10% 3p21. 31 (46,827,053- 




del 13q 14.3 






50% 3p21.1(52,038,099-52,553,933)x1 b , >90% 


50,133,778)x1 b , 10% 3p21 .1 (52,038,099- 










1 1q14.3q24.3(91, 864,846-1 27,805,084)x1, >90% 


52,553,933)x1 b , 80-90% dell 1 q14.3q24.3 










1 3q1 4.1 1 (40,389,1 70-40,761 ,341 )x1 b , > 90% 


(91, 864,846-1 27,805,084)x1, 90% 13q14.11 










1 3q 1 4.1 1 (41 ,687,91 2-41 ,826,01 6)x1 b , > 90% 


(40,389,1 70-40,761,341 )x1 b , 90% 13q14.11 










1 3q 1 4.2q 1 4.3(47,5 1 5,748-5 1 ,333,756)x1 


(41,687,91 2-41 , 826,01 6)x1 b , 90% del13q14.2q14.3 












(47,51 5,748-51, 333,756)x1 


CLL108 


None 


Un mutated 


1 


46,XX, arr 3q22.3(1 37,678,290-1 37,832,637)x1 b , 


46,XX, arr 70% 6q1 4.1 q22.33(81, 463,977- 










80% 6q 1 4.1 q22.33(81 ,463,977-1 29,05 1 ,286)x1 , 


1 29,051, 286)x1, 70% 12(0-1 32,1 46,663)x2 hmz, 










90% 1 1q12.3qter(62,398,575-134,445,626)x2 


30% 17pterp1 3.1 (0-8,033,41 9)x1, 12p12.2q21.31 










hm7 RO% 1 9ntprn1 1 HCl-lCl hm7 Rfi% 

i ii i iz., o\j /u i z.|jlci yj i z. .z. \\j zu,j JU/T^J/ y AZ. 1 II 1 IZ., ou /u 


Dfi 4S7-R4 9S1 9RAW3 nn Innnpr Hptprtahlp 










1 2p1 2.2q21 .31 (20,356,457-84,25 1,286)x3, 80% 












1 2q21 .31 qter(84,251 ,286-1 32,1 46,663)x2 hmz 




CLL003 


de11q22.3 


Un mutated 


2 


46,XY, arr 1 1q22.3(1 04,1 58,306-1 09,733,1 34)x1, 


46,XY, arr 5q12.1(58,925,836-58,976,604)x1 




del 13q 14.3 






11q23.1q23.2(1 10,556,535- 












1 1 4,089,91 0)x1 ,1 1 q23.3(chr1 1:11 5,527,657- 












1 16,182,933)x1, 13q 14.3(49,038,530- 
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{Continued) 










Patient 


Known 


IGHV gene 


CRS 


CNAs and cnLOH regions in pre-treatment sample 


Differences at relapse 


no. 


CNAs 


status 
















50,906,239)x1 , 1 7q21 .33q22(45,720,750- 
49,01 7,380)x2 hmz c 














CLL006 


tri12 


Un mutated 


2 


47,XY, arr 90% 4q31. 1(1 53,631 ,568- 
153,935,927)x1, 70% 12(0-1 32,1 46,663)x3, 
1 7p1 2(1 3,1 1 8,349-1 3,1 71 ,201 )x3 b 


None noted 


CLL044 


del13q14.3 
de17p13.1 


Mutated 


2 


46,XY, arr 46,XY, arr 2q36.2 (225,094,102- 
225,61 7,560)x1, -20% 9p23p2 1.1 (1 1, 986,1 63- 
32,324,220)x1 b , 50% 9q13qter(70,223,358- 
1 39,946,600)x3, 1 3q 1 4.3(49,5 1 5,641 - 
51 ,327,1 98)x0, 1 7pterp1 1 .2(0-21 ,430,683)x1 , 
1 9q1 3.1 2q1 3.2(41 ,320,887-43,41 2,844)x2 hmz c 


46,XY, arr -30% 9p23p21. 1(1 1,986,1 63- 
32,324,220)x1 b , <20% 9q13qter(70,223,358- 
139,946,600)x3 b 


CLL056 


del13q14.3 


NK 


2 


46,XY, arr 80% 5p1 3.3(37,408,772-37,784,1 35)x3 b , 

80% 9pter9p11.2(0-43,254,687)x2 hmz c , 70% 

1 3q 1 4.1 2q 1 4.2(45,377,644-45,760,335)x3, 70% 

1 3q 1 4.2(46,403,477-47,574,445)x3, > 90% 

1 3q 1 4.2(47,583,91 7-48,099,556)x1 , 90% 

1 3q 1 4.2q 1 4.3(48,7 1 9,245-50,523,975)x1 , 70% 

1 3q 1 4.3qter(50,523,975-1 1 4,1 21 ,252)x3 


46,XY, arr 30% 4pterp1 5.1 (0-35,1 66,759)x1 b , 90% 
9pter9p11.2(0-43,254,687)x2 hmz c , -50% 
Xq26.2q28 (1 31 ,629,636-1 54,058,61 7)x2 b 


CLL063 


de11q22.3 
del13q14.3 


Un mutated 


2 


46, XX, arr 1 1q1 4.1 q24.1 (77,652,729- 
1 22,488,1 38)x1 , 1 3q 1 4.1 1 (39,758,892- 
39,909,885)x0 b , 1 3q1 4.3(49,425,974- 


None noted 










50,389,21 0)x0 b , 1 3qcenqter(1 9,075,593- 
114,142,980)x2 hmz c 














CLL079 


None 


Un mutated 


2 


None noted 


46, XX, arr 6p25.3(746,61 8-789,781 )x3 b , 
6p22.3(1 9,723,828-1 9,765,1 53)x3 b , - 40% 
6q 1 6.3(1 01 ,257,1 01-101 ,307,223)x3 b , 
6q21 (1 09,405,342-1 09,444,628)x1 b 


CLL082 


del13q14.3 


Mutated 


2 


46, XX, arr 13q 13.2-21 .2(33,726,200- 

53,61 9,700)x1 , 1 3q1 4.3(49,461 ,448-50,635,061 )x0 


None noted 


CLL084 


None 


Un mutated 


2 


46,XY, arr 4q34.1 (1 72,446,71 7-1 72,482,835)x1 b 


46,XY, arr < 10% 6pterp22.3(0- 
17,028,648)x1 b ,<10% 10q24.31qter (102,668,100- 
135,374,737)x1 b , <10% 11q14.1q23.2 (78,372,164- 
1 1 4,844,825)x1 b , < 1 0% 1 7pterp1 1 .2(0- 
1 8,8771 65)x1 b , <10% 18pterpcen (0- 
13,907,055)x1 b , <10% 1 8q21.1q21. 2(41, 130,378- 
50,155,862)x1 b , <10% 18q21 .33q22.1 (59,798,490- 
62,065,41 8)x1 b , <10% dell 8q22.2qter 
(66,972,378-76,11 7,1 53)x1 b 


CLL086 


None 


Mutated 


2 


46,XY, arr 3p26.3p26.1 (0-7,571, 446)x1, >90% 
Xq27.1 qter(1 39,443,330-1 54,582,606)x2 


46,XY, arr 11q14.3 (89,597,760-89,690,254)x3 b , 
1 2q24.23(1 1 8,81 1 ,967-1 1 8,875,247)x3 b , 
22q1 2.1 (26,427,496-26,490,006)x3 b 


CLL087 


tri12 


Un mutated 


2 


47,XY, arr 80% 12(0-1 32,1 46,663)x3, 
Xpl 1 .3p1 1 .4(41 ,074,455-43,983,367)x0 


None noted 


CLL088 


del 13q 14.3 


Mutated 


2 


46,XX, arr >90% 1 3q 1 4.1 2q 14.3(45,81 7, 
763-49,547,499)x1 


None noted 


CLL095 


tri12 


NK 


2 


47,XY, arr 80% 2q14.3qter(1 29,537,1 70- 
242,497,808)x2 hmz, 6q 1 4.1 q21 (79,968,009- 
1 06,91 7,920)x1, 70% 12(0-1 32,1 46,663)x3, >90% 
20q11.2qter (30,155,027-62,435,964)x2 hmz c , 
22q1 1 .23(35,01 0,463-35,020,688)x0 b 


46,XY, arr 2p14(chr2:67,974,461-68,005,924)x3 b 


CLL102 


tri12 


Un mutated 


2 


47,XY, arr 80% 12(0-1 32,1 46,663)x3 


None noted 


CLL103 


tri12 


Un mutated 


2 


47,XX, arr 1q25.3(1 83,360,633-1 83,397,1 69)x3 b , 
> 90% 2q22.1 q31 .1 (1 41 ,094,077-1 71 ,777,386)x2 
hmz, >90% 9p24.1p22.3(9,000,336-14,313,346)x2 
hmz, 10q26.3(1 35,089,371 -135,374,737)x3 b , 80% 
12(0-1 32,1 46,663)x3, >90% 
1 3q33.3q34(1 08,839,643-1 1 2,01 4,298)x2 hmz, 
70% 22q1 1 .21 q1 2(1 7,295,655-1 9,791, 274)x3 


None noted 


CLL106 


del17p13.1 


Un mutated 


2 


46,XY, arr 90% 1 p34.2(41, 538,823-42,1 12,41 0)x1, 
90% 6p25.3p25.1 (1,376,1 35-6,1 73,21 4)x1, 90% 
6p25.1 p24.3(6,471,929-7,578,657)x1, 90% 
6p24.3p24.1 (10,1 55,458-1 1,223,334)x1, 90% 
6p24.1 p22.3(1 2,754,522-1 8,670,995)x1 , 90% 
6p21 .2(38,889,727-39,544,026)x1, 90% 
6p21 .1 p1 2.3(43,689,563-46,299,596)x1,90% 
7p1 5.3(23,346,71 0-24,400,934)x1, 90% 9q33.3 
(127,381, 099-1 27,900,525)x1, 90% 12p13.33 
(977,044-1 ,959,491 )x1 , 1 7pterpcen(0- 
22,1 25,1 62)x1 , 1 8pterp1 1 .31 (0-6,572,1 04)x1 , 
1 8p1 1 .31 pi 1 .23(6,802,509-7,333,763)x3, 


46,XY, arr 80% 1 p34.2(41, 538,823-42,1 12,41 0)x1 , 
60% 3p26.3p26.1(0-8,686,575)x1, 70% 
3p25.1 p22.3(1 2,21 4,891 -34,481 ,950)x1 , 70% 
3p21 .33p21 .31 (43,708,1 88-46,679,233)x1 , 70% 
3p21 .1 p1 4.2(52,51 5,81 3-61 ,534,81 1 )x1 , 70% 
3p1 4.1 p1 3(68,801 ,809-74,909,820)x1 , 60% 
8pterq11.21(0-49,228,092)x1, 60% 8q12.1q13.2 
(57,01 1,940-70,263,279)x1, 60% 8q23.1 
(1 06,323,747-1 06,959,523)x1, 90% 1 7pterpcen(0- 
22,1 25,1 62)x1, 90% 18pter18p1 1.31(0- 
6,572,104)x1, 90% 18p1 1.31 p1 1.23(6,802,509- 
7,333,763)x3, 90% 1 8p1 1 .22p1 1 .21 (8,926,107- 
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Table 3 {Continued) 



Patient Known 
no. CNAs 



IGHV gene 
status 



CRS CNAs and cnLOH regions in pre-treatment sample Differences at relapse 



CLL145 de17p13.1 Unmutated 2 



1 8p1 1 .22p1 1 .21 (8,926,1 07-1 3,281 ,880)x1 , 
1 8p1 1 .21 (1 3,536,937-1 5,348,275)x1 , 90% 
22q1 1 .21 (1 6,71 2,51 0-1 7,055,91 4)x1 b , 90% 
22q1 1.23q12.1(23,409,925-27,331,789)x1, 90% 
22q1 2.1 ql 2.2(27,573,41 0-27,91 8,1 40)x1 b 
46,XY, arr <10% 3p25.1 pi 4.3(1 2,51 6,876- 
55.709725)x1 b , <10% 5q14.3q23.3(87,285,075- 
129,855,075)x1 b , >90% 6p24.3p22.3(1 0,220,342- 
1 6,992,531 )x2 hmz b , 6q1 6.3(1 00,876,399- 
101,445,547)x1, 80% 8q22.3q24.21 (106.058,337- 
1 28,286,81 8)x2 hmz, 80% 8q24.21 (128,286,81 9- 

1 28,297,901 )x3 b , 80% 8q24.21-q24.3(1 28,297,902- 17p13.3p1 1.2(0-1 8,485,225)x2 hmz b , -20% 



1 3,281, 880)x1, 90% 18p1 1.21 (13,536,937- 
15,348,275)x1, 80% 22q1 1.21 (16,71 2,510- 
1 7,055,91 4)x1 b , 80% 22q1 1.23q 12.1 (23,409,925- 
27,33 1,789)x1, 80% 22q 12.1 q1 2.2(27,573,410- 
27,91 8,1 40)x1 b 

10-20% 3p25.1p1 4.3(1 2,51 6,876- 55.709725)x1 b , 
10-20% 5q14.3q23.3(87,285,075- 
129,855,075)x1 b , 50% 8q22.3q24.21 (106.058,337- 
1 28,286,81 8)x2 hmz, 50% 8q24.21 (128,286,81 9- 
1 28,297,901 )x3 b , 50% 8q24.21-q24.3 
(1 28,297,902-1 47,274,826)x2 hmz b , <10%-40% 



12q14.1q14.2(56,925,673- 
16q22.3q23.1 (72,290,299- 
<10%-50% 17p13.3p1 1.2(0- 
- 30% 1 8q 1 1 q23(22,045,087- 
80% 22q1 1.23(,23,991,557- 



147,274,826)x2 hmz b 
62,754,801 )x2 hmz b , 
76,525,557)x2 hmz b , 
18,485,225)x2 hmz b , 
76,11 5,1 72)x2 hmz b , 
24,250,622)x1 b 

CLL156 None NK 2 46,XY, arr 80% 5q 23.3 (1 24,282,026- 

124,531,381)x3 b , -10% 

1 9q1 3.1 2q1 3.1 3(41 ,440,201 -43,1 49,438)x1 hmz b 
ARR003 del13q14.3 Unmutated 2 46,XY, 13q14.2q 14.3(49,983,042-51, 33 1,787)x1 b 

CLL071 del13q14.3 Mutated 3 46,XY, arr < 10% 13q13.3q34(36,483,665- 



CLL081 del 13q 14.3 
de17p13.1 



Unmutated 3 



CLL092 de11q22.3 
del 13q 14.3 



Unmutated 3 



CLL094 de17p13.1 Unmutated 3 



1 1 4,1 23,1 25)x1, 70% 1 3q14.2q 14.3(48,7 13,406- 
50,479,598)x1 b , 1 3q1 4.3(49,229,284-50,423,376)x0, 
30-70% 17pterp11. 2(0-21, 100,678)x2 hmz 
46,XY, arr <10% 2pterp1 4(0-67,98 1,586)x3 b , 90% 
2q 1 1 .1 q 1 1 .2(95,597,687-97,676,31 4)x1 , 90% 
2q 1 4.2q 1 4.3(1 21 ,572,920-1 26,53 1 ,524)x1 , 90% 
2q 1 4.3(1 27,036,406-1 28,569,25 1 )x1 , 90% 
2q21 .2q24.1 (1 34,752,1 38-1 55,542,922)x1 , 40% 
8q22.2qter(1 00,925,1 55-1 46,263,890)x3, 90% 
9pterp1 1.2(0-33,31 8,81 2)x1, 90% 9p13.3p11.2 
(33,83 1,777-43,548,708)x1, >90% 9q21.13 
(74,372,294-74,5 16,002)x1 b , <10% 10q23.2q23.31 
(89,306,467-93,720,035)x1 b , 90% 11q12.2q13.1 
(61,231,809-64,790,450)x1, 70% 1 2q1 3.1 1 ql 3.1 2 
(47,308,669-48,029,1 06)x3, 50% 13q11q13.1 
(17,920,393-32,659,455)x1, 50% 13q14.3 
(48,898,975-50,805,21 9)x1, 90% 1 7pterpcen(0- 
22,1 25,1 62)x1, 70% 1 7qcenqter(22,71 3,949- 
78,637,1 98)x3; 90% 19pterp1 3.2(0-1 0,071, 486)x1; 
70% 1 9p1 3.2p1 3.1 3(1 0,071 ,486-1 3,41 9,854)x3 

46,XY, arr 2p16.1(55,330,268-55,352,382)x1 b , 
1 1 q22.3q24.3(1 03,008,248-1 27,852,266)x1 , 
1 3q1 4.3(49,500,283-50,255,794)x1 , - 20% 
20q1 1.22qter(32,772,470-62,385,675)x2 hmz c 
46,XY, arr 30% 2q33.2q36.3(206,690,791- 
230,780,953)x1, 40% 17pterp1 1.2(0-20,067,160) 
x1 b , 60% 17p1 1.2(20,067,1 61 -21, 471, 049)x1 b 



1 8q1 1 q23(22,045,087-76,1 1 5,1 72)x2 hmz c , 



None noted 



46,XY, arr 50% 1 1q13.1q25(66,992,701- 

134,177,292)x2 hmz b , 90% 

1 3q 1 4.2q 1 4.3(49,970,598-5 1 ,333,659)x1 fc 

46,XY, -20% 1 3q13.3q1 4.2(36,446,865- 

114,123,122)x1 b . 



46,XY, arr -50% 2pterp1 4(0-67,98 1,586)x3 b , 
>90% 2q1 1.1 q1 1. 2(95,597,687-97,676,31 4)x1, 

> 90% 2q1 4.2q1 4.3(1 21 ,572,920-1 26,531 ,524)x1 , 

> 90% 2q1 4.3(1 27,036,406-1 28,569,251 )x1 , 

> 90% 2q21 .2q24.1 (1 34,752,1 38-1 55,542,922)x1 , 
- 1 0% 2q31 .1 (1 74,294,1 44-1 74,955,1 41 )x1 b , 

- 30% 2q33.1 q37.1 (204,070,804-231 ,244,281 )x1 b , 
90% 8q22.2qter(1 00,925,1 55-1 46,263,890)x3, 
>90% 9pterp1 1.2(0-33,31 8,81 2)x1, >90% 
9p13.3p1 1. 2(33,83 1,777-43,548,708)x1, 40% 
1 0q23.2q23.3 1 (87,41 9,694-1 07,982,339)x1 , 
>90% 11q12.2q13.1(61,231,809- 
64,790,450)x1 , - 1 0% 1 3q 1 1 q 1 3.1 (1 7,920,393- 
32,659,455)x1, <10% 1 3q1 4.3(48,898,975- 
50,805,21 9)x1 b , -30% 1 6p1 3.3(3,407,330- 
4,492,425)x1 b , >90% 1 7pterpcen(0- 
22,1 25,1 62)x1 , > 90%1 9pterp1 3.2(0- 
1 0,071, 486)x1 

46,XY, -10% 1 1q22.1 (101, 354,200- 

101,514,796)x1 b , arr 40% 

20q1 1.22qter(32,772,470-62,385,675)x2 hmz 

46,XY, arr 20% 2pterp21 (0-45,928,1 91 )x1, 20% 
del2p1 6.32p1 6.1 (48,61 7,246-59,276,760)x1 , 20% 
2p1 3.2p1 1 .2(73,382,585-83,989,328)x1 b , 
2q33.2q36.3(206,690,791 -230,780,953)x1 (no 
longer detectable); 20% 6p24.3p25.2(3,1 38,439- 
9,391,1 13)x1, 20% 6p22.3(1 6,521,1 54- 
19,237,784)x1, 20% 6p21. 31 (35,1 54,526- 
36,820,525)x1 b , 20% 6p21 .2(3 7,804,382- 
38,908,035)x1 b , 20% 6p21 .1 q1 2(48,685,550- 
67,685,569)x1 b , 20% 6q1 3q1 5(74,279,874- 
88,065,764)x1, 20% 6q21q22.1 (1 1 1, 800,880- 
116,278,896)x1 b , 20% 6q22.33(1 27,739,685- 
130,160,593)x1, 20% 7q21 .1 1 (80,379,709- 
85,591,71 1)x1, 20% 7q31.1q32.1(1 1 1, 820,986- 
127,193,285)x1, 20% 7q32.3qter(1 30,945,628- 
157,796,605)x1, 20% 8pterp12(0-36,598,375)x1 b , 
20% 1 0q24.1 qter(98,893,384-1 35,085,880)x1 , 
50% 11q13.1q25 intermittent(63,515141- 
131,026,292)x3 b , 70% 17q1 1.2(24,979,666- 
25,506,441 )x1, 70% 1 7q1 1 .2(27,442,088- 
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Table 3 {Continued) 



Patient Known 
no. CNAs 



IGHV gene 
status 



CR5 CNAs and cnLOH regions in pre-treatment sample Differences at relapse 



CLL112 del13q14.3 NK 



46,XYarr 90% 2q35(218,476,667-220,204,673)x1, 
90% 2q36.3q37.1 (230,239,060-231 ,446,51 3)x1 , 
8q21 .1 1 q21 .1 3(78,1 81 ,345-80,471 ,1 88)x2 hmz c 
, > 90% 1 3q 1 1 .2qter(20,381 ,113-11 3,908,097)x2 
hmz, 13q1 4.3(49,494,577-50,238,31 6)x0, 80% 
1 6p1 2.1 (21 ,860,290-22,350,261 )x1 



27,692,902)x1 b , 60% 1 9p13.3p1 3.1 1(5,676,493- 
18,095,588)x3, 70% 1 9p1 3.1 1 p1 2(1 8,201,910- 
24,337,497)x1 

46,XYarr 90% 2q35(218,476,667-220,204,673)x1 no 
longer detectable, 90% 2q36.3q37.1 
(230,239,060-231,446,51 3)x1 no longer 
detectable, >90% 13q1 1.2qter(20,381,1 13- 
1 1 3,908,097)x2 hmz no longer detectable, 
13q14.3(49,494,577-50,238,316)x0 no longer 
detectable 



Abbreviations: CNA, copy number alteration; cnLOH, copy neutral loss of heterozygosity; CRS, clinical risk score; hmz, homozygous; IGVH, immunoglobulin 
variable heavy chain gene; NK, not known. All co-ordinates are hg36 and determined using OncoSNP except '° (see above). 'x2 hmz' indicates cnLOH. a Using 
size thresholds of 20 kb for CNAs (unless in a known region of interest) and 2 Mb for cnLOH. Percentages represent levels of mosaicism estimated from B-allele 
frequency values. Where a percentage is not noted, then all cells appear to carry the anomaly. b Co-ordinates determined using Nexus (because of variable CNA 
boundaries less well defined by OncoSNP or <500kb in size or below the level of detection of OncoSNP). c Called by OncoSNP as possible germline variant 
and manually verified as probable somatic event. Co-ordinates determined using Nexus. 
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Figure 2. OncoSNP output showing clonal expansion events involving the chromosome 13q14.3 MDR and percentages of cells involved in 
pre-treatment and relapse samples from CLL080 and CLL096. The red lines indicate deletion events, whereas the magenta lines indicate allelic 
imbalance/cnLOH. In CLL080, 30% cells carry a large 13q14.1 1q14.3 deletion and an additional ~20% cells carry the smaller MDR deletion at 
pre-treatment. At relapse, ~90% cells carry the larger deletion. For CLL096, there is a homozygous deletion of 13q14.3, and cnLOH of the 
entire chromosome in ~40% cells at pre-treatment. This cnLOH expands to ~90% cells at relapse. 



involved in cancer proliferation and survival. A fourth region affects DNA repair genes (BTBD12, DNASE1), TRAP1, a mitochondrial 
NFIB involved in the NFKB pathway and two genes of no obvious chaperone and regulator of apoptosis, as well as CREBBP, a 
CLL related annotation (TYRP1, MPDZ). A further region included frequently mutated gene in B-cell lymphomas. 
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Figure 3. OncoSNP output showing clonal evolution in pre-treatment and relapse samples from CLL081 . This patient shows a reduction in the 
percentage of cells carrying the 13q abnormalities whilst other CNAs expanded (8q22.2qter and 10q23.2q23.3) or were newly identified at 
relapse (2q33.1q37.1 and 16p13.3). 



When we performed pathway analysis of the 546 genes in 
the MDRs/MORs affected by recurrent and emerging/expanding 
CNAs, we identified 17 over-represented, statistically significant 
and independent pathways of which eight contained 55 cancer- 
related genes (phagosome P-value = 0.0024; apoptosis P-value = 
0.0079; small-cell lung cancer P-value = 0.0232; prostate cancer 
P-value= 0.0284; pancreatic cancer P-value = 0.0355; chronic 
myeloid leukemia P-value = 0.0414; wnt signaling P-value = 0.0436; 
cell cycle P-value = 0.0497) (see Supplementary Table 10, online). 

DISCUSSION 

Our study is the first genome-wide array based analysis of CNAs/ 
cnLOHs that both characterizes and quantifies the proportion of 
subclones carrying genomic changes before treatment and at 
subsequent relapse in cancer. The results demonstrate that 
(i) many CNAs/cnLOHs in CLL are recurrent and therefore non- 
random events that expand over time owing to Darwinian 
selective pressure, (ii) selected genes identified in both well- 
recognized and newly defined MDRs/MORs present plausible 
candidates for driving disease progression and include those 
involved in mature B-cell development, DNA damage response 
(DNA repair, apoptosis and autophagy), tumor progression 



and familial B-CLL; pathway analysis was consistent with this. 
The identification of genes in this way will be informative for 
focused sequencing strategies, (iii) In our cohort there is strong 
supportive evidence that genomic complexity is associated with 
poor risk disease. This adds strength to previous studies that 
indicated a link between genomic complexity or clonal evolution 
to survival in CLL 30,36,38 (iv) OncoSNP is a powerful computational 
statistical tool that alters fundamentally our ability to investigate 
clonal architectures and tease apart complex clonal dynamics. 
Furthermore, our results demonstrate the importance of devel- 
oping comparable algorithms for application to next generation 
sequencing data. Although it is likely that in the longer term array 
technology will be replaced by whole genome sequencing, this 
study contributes to our understanding of the significance of 
CNAs/cnLOH in cancer progression. Finally, our results do not 
support the existence of a single leukemia propagating cell 
population but instead, shed light on the considerable genomic 
heterogeneity of cells driving disease progression. 

Clinically, the existence of multiple, genetically distinct, 
subpopulations that escape therapeutic intervention presents 
formidable challenges for the development of effective treat- 
ments for patients with relapsed refractory B-CLL. Therefore, 
detailed characterization of the molecular basis of the condition 
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Figure 4. Schematic representation of the possible clonal architecture pre-treatment and at relapse (pale blue and mid-blue shaded areas, 
respectively) for patients (a) CLL106 and (b) CLL081. Proportions below 10% and above 90% and cell populations without identifiable CNAs/ 
cnLOHs are inferred. The precision of the percentage contribution of subclones with CNAs/cnLOH is ±10%. Note that for CLL081, the 
proportion of subclones with the 13q14.3 loss decreases at relapse and the 1 2q 1 3.1 1q13.12 and 1 9p1 3.1 3p1 3.2 gains are no longer observed 
in one subclone. 



and elucidation of the underlying mechanisms of clonal diversity 
will be essential for future targeted selection of effective 
therapeutic agents. 
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